Disorder effects in the quantum Heisenberg model: An Extended Dynamical 

mean-field theory analysis 
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We investigate a quantum Heisenberg model with both antiferromagnetic and disordered nearest- 
neighbor couplings. We use an extended dynamical mean-field approach, which reduces the lattice 
problem to a self-consistent local impurity problem that we solve by using a quantum Monte Carlo 
algorithm. We consider both two- and three-dimensional antiferromagnetic spin fluctuations and 
systematically analyze the effect of disorder. We find that in three dimensions for any small amount 
of disorder a spin-glass phase is realized. In two dimensions, while clean systems display the proper- 
ties of a highly correlated spin-liquid (where the local spin susceptibility has a non-integer power-low 
frequency and/or temperature dependence), in the present case this behavior is more elusive unless 
disorder is very small. This is because the spin-glass transition temperature leaves only an interme- 
diate temperature regime where the system can display the spin-liquid behavior, which turns out to 
be more apparent in the static than in the dynamical susceptibility. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.50.Ee, 75.40.Gb 



I. INTRODUCTION 



Quantum magnetism is an important topic of mod- 
ern solid state physics: Not only do magnetic phases 
take place in the phase diagram of many correlated 
electron systems, like high-temperature superconducting 
cupratesi and heavy fermions^, but many different ma- 
terials display magnetic properties, which change upon 
varying parameters such as doping, temperature, disor- 
der, or pressure. A few examples may (very partially) 
illustrate the enormous variety of behaviors that one 
encounters in this field: Upon increasing doping the 
lightly doped La2_ x Sr x Cu04 passes from an antiferro- 
magnetic (AF) to a spin-glass (SG) phased; the intro- 
duction of non-magnetic impurities in quantum paramag- 
nets like SrCu203, KCuCi3, CuGe03, gives rise to local- 
moment formation 3 -; the glassy phase of the disordered 
dipolar-coupled magnet LiHo x Yi_ x F4 may be driven to 
other magnetic phases by means of a magnetic field 4 ; 
SrCr9 P Gai2-9pOi9 is a S — 3/2 magnet with a Kagome 
lattice structure, which above the (low) SG freezing tem- 
perature displays the spin-liquid behavior 5 '- typical of 
strongly correlated spin systems^. From this large va- 
riety we extract three essential ingredients, which are 
commonly and physically relevant in quantum magnetic 
materials: disorder, frustration, and fluctuations (both 
thermal and quantum) . The separate role of these physi- 
cal mechanisms and their interplay in determining differ- 
ent physical properties is obviously a relevant issue in the 
field of quantum magnetism. Since one is facing a huge 
variety of systems with different structures, properties 
and phase diagrams, the search of common fundamental 
mechanisms mostly relies on the analysis of schematic 
models to highlight the deep physical effects. For in- 
stance some work has been devoted to the interplay be- 



tween disorder and fluctuations^^, analyzing the role 
of quantum spin fluctuations in driving a spin-glass (SG) 
phase into a paramagnetic (or spin-liquid) phase upon 
varying spin size and/or temperature. Within the same 
framework we consider in this paper a "simple" quantum 
AF Heisenberg model for S = 1/2 spins and we intro- 
duce a random magnetic coupling. By solving this Disor- 
dered Quantum Antiferromagnetic Heisenberg (DQAFH) 
model within the Extended Dynamical Mean-Field The- 
ory (EDMFT) 11 ! 12 ' 13 , we aim to clarify the interplay be- 
tween AF and random magnetic couplings as well as the 
role of dimensionality in the instability from a paramag- 
netic to a SG phase and the relevance of quantum fluctu- 
ations in suppressing the critical temperature. Through- 
out the present work we will always remain in the para- 
magnetic phase with neither AF- nor replica-symmetry 
breaking and we will determine the conditions for these 
(second-order) instabilities to take place. 

The EDMFT techniqueii^ii 3 - is an extension of the 
standard Dynamical Mean-Field Theory^, where the lo- 
cal (spin or charge) degrees of freedom are coupled to a 
self-consistent bosonic bath. In our case the Heisenberg 
model is mapped onto a single quantum spin embedded 
in a bosonic bath which represents the surrounding spin 
excitations as a fluctuating magnetic field. This approach 
has some definite limitations: It is not suited to keep full 
account of the spin character of the surrounding degrees 
of freedom (for instance it misses the possibility of sin- 
glet formation between the single spin and the surround- 
ing ones) and it neglects the momentum dependence of 
the spin self-energy thereby preventing the occurrence of 
spatial anomalous dimensions. This latter limitation is 
more severe at low temperatures, where spatial correla- 
tions between the spins become relevant. As also dis- 
cussed in Section VI the consequence of this is the ap- 
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pearance of well-know n 13 ' 15 low-temperature instabilities 
occurring as an unphysical sign in the spin self-energy. 
Also related to this poor treatment of spatial correla- 
tions (and the related lack of anomalous dimensions) is 
the appearance of spurious first-order transitions, where 
second-order are instead expected^ 6 -. Nevertheless, it is 
commonly recognized 16 that the EDMFT approach is a 
valuable tool to investigate the dynamics of spin systems 
and it makes an important first step toward the inclusion 
of spatial correlations between the spins. 

Our analysis of the DQAFH not only has a broad valid- 
ity and addresses general theoretical issues, but it is also 
of pertinence for those numerous AF materials, where 
disorder can induce competition between antiferromag- 
netism (possibly having a correlated anomalous charac- 
ter) and a spin glass (SG) phase, or can induce local 
moment formation with low-temperature Curie-like de- 
pendencies. The role of low dimensionality in determin- 
ing the relevance of magnetic fluctuations is also investi- 
gated. 

The paper is organized as follows. In Section II we 
present the model and the EDMFT technique. Section 
III presents the criteria necessary to detect the establish- 
ment of broken-symmetry phases or the spin-liquid be- 
havior. The numerical results for static properties in the 
three- and two-dimensional case are reported in Section 

IV, while the dynamical behavior is discussed in Sect. 

V. In Section VI we discuss our results and present our 
conclusions. 



II. THE MODEL AND THE TECHNIQUE 

A. The model 

We consider the Disordered Quantum Antiferromag- 
netic Heisenberg (DQAFH) model, given by the following 
Hamiltonian: 



i<j 



3 • 



(1) 



AF = JiA F /\fz = J/Vz [with Fourier 



Here Si is a spin— 1/2 operator at the i— th site of a lat- 
tice of coordination z. The nearest-neighbor magnetic 
couplings Jij = jf- F + J® are the sum of a translation- 
ally invariant part 

transforms </(q)] and a disordered part = J-P/yfz. 
The translationally invariant couplings are related to 
the spin-wave spectral density p(e) = J2q^( e ~ ^(q))- 
The disordered parts J® are random couplings with a 
Gaussian probability distribution P(J®), with zero mean 

value (Jij)dis = and {J-jJ^dis = bi,i&j,k J D [(••*)<«« 
represents the averaging over the disorder distribution] . 

In the following, we will perform a large- z expansion, in 
order to study the spin-wave correlations within a single- 
site effective model. This derivation has been described 
in Ref. for a clean Heisenberg model. Here, the pres- 
ence of disorder requires a specific preliminary treatment. 



1. Disorder and replicas 

The free energy is averaged with respect to the prob- 
ability P(jfj). We use the well-known replica trickii, 
which relies on the fact that the free energy T can be 
calculated from the partition function Z using the rela- 
tion 

f3T = (\ogZ) dls = lim - \og{Z n ) dls . (2) 

In order to calculate {Z n ),i is we introduce n replicas of 
the system, each one being represented by an index a 



(z n ) 



dis 



n 



+ OC 



dJ°P(jP)Tr 



exp 



a=l 



where H a is the disordered Heisenberg Hamiltonian Eq. 
([1]) acting on the replica a. In Eq.Q, the trace is applied 
to all the operators Sf of the replicated Hilbert space (see 
Appendix A). 



2. Extended Dynamical Mean Field 

In the standard DMFT technique^, the electronic de- 
grees of freedom of a specific single site are coupled to 
a fermionic effective bath. On the other hand, in the 
EDMFT spin (or charge) degrees of freedom are coupled 
to an external bosonic bath simulating the spin/charge 
degrees of freedom of the rest of the lattic o 11 ' 12 ' 13 . The 
condition that the single specific site must be equivalent 
to all other sites of the lattice is then implemented by 
means of selfconsistency equations. We are now going to 
suitably adapt this procedure to the disordered spin sys- 
tem represented in Eq. ([1]). The first step is standard: 
We single out a specific site of the system (arbitrarily 
chosen as the origin); the disorder (J 9 couplings) will 
be formally averaged and the trace will be taken with 
respect to the operators Sf with i ^ 0. For any value of 
z, one obtains within this procedure an effective local ac- 
tion for the local spins Sg . The local action is calculated 
with a cumulant expansion 18 : In the large-z limit, all the 
cumulants vanish, except the first two terms. One de- 
pends linearly on Sf and corresponds to the static mean 
field. The second is quadratic and reflects the dynamical 
fluctuations. Here, we study only the paramagnetic and 
replica-symmetric solution with (Sf) = 0: The transition 
to an AF phase or to a SG phase will only be analyzed 
from the paramagnetic phase, considering the criteria for 
a second-order transition (see Sect. III). 

The local effective action is thus cast into the following 
form: 



A=- 



~ f dr f dT'Y J ^W{T')K, aa, {T- 
l Jo Jo 



r'), 
(4) 
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with 

JC aa '(T-r') = l]T(J 0l J j0 (S«(T)Sf {r')) camtv ) dls ■ 

ij 

In this expression for the Kernel, as in the following ones, 
the summation with respect to i and j concerns the near- 
est neighbors of the site 0, while a, a' are replica indexes. 
The cavity correlation functions (• • • Y amty correspond to 
the correlation functions calculated with a cavity instead 
of the site (see Appendices A and B). The coupling fac- 
tors JoiJjo in Eq. ([5]) precisely invoke the site (of the 
cavity) . 

Using the arguments of Appendix B, the averaging over 
the disorder can be factorized 

(j 0l j j0 (^(r)sf(T')r mt y) dls ^ 

(JoiJ^disiiSt^sfiT'))^}^ . (6) 

Since we are concerned with (the instabilities of) the 
paramagnetic phase we only consider averages (with and 
without cavities) leading to diagonal quantities both in 
the spin components and in the replica indexes. All the 
following arguments could be carried out keeping the ma- 
trix character of the spin correlation functions, but we 
adopt this simplified presentation to put more clearly in 
evidence the role of disorder without unnecessary formal 
complications. In this way all thermally averaged quan- 
tities are scalar quantities and accordingly 



JC(t-t') = /C QQ (r 



(7) 



and one can introduce the (diagonal, i.e., scalar) spin 
susceptibility 

Xij (r - /) = (S?(t)S?(tO>/3 = (Sr(r)S^(r')) . (8) 

\(i is a spin-component index] . We also define the disorder 
averaged susceptibilities 



X, 



(x 



ij /dis 



(9) 



with a similar definition for the average cavity suscepti- 
bilities Xij VltV ■ Now, invoking Eq. ((6| in Eq. ((5|), fixing 
the unnecessary replica index and dropping it, we find 
for the paramagnetic phase 



£( T ) = - ^2{JoiJoo)disXt] mtV {T) 



(10) 



In the large— z limit, the averaged cavity susceptibility 
can be expressed in terms of the full physical (i.e. without 
cavity) averaged susceptibility (see Appendix B) 

%j mtV (Un) = ^(Un) -XrtMXojM/XooM ■ (11) 

where u> n is a Matsubara frequency. The Kernel (fTU)) is 
then given by the self-consistent relation 

£ = - ( J (K Fj jo F + J D^j) [Xij ~ XioXoi/Xoo] . 

ij 

(12) 



where the explicit Matsubara frequency dependency 
has been dropped for clarity. Using standard 
consideration o 11 ' 14 ' 19 of power counting estimates of the 
correlation function dependence on l/z one can cast the 
part proportional to J^S.^ in a fully local form 



K 



z ^ 

K-AF 



JtM F Jjo F (Xij XioXoj/Xoo) + J 



r>Xi 



(13) 



Although the disorder obviously enters the averaged x's, 
it is important to notice that in Eq. (fT3"|) the kernel splits 
in a part only involving ordered AF couplings and the 
generic \ij ~ XioXoj/Xoo combination, and a term pro- 
portional to J|> and to the purely local susceptibility Xu ■ 
This decomposition of the kernel into an AF term involv- 
ing non-local correlations and a D term with only local 
correlations is one of the noticeable results of this work. 
Therefore, as far as disorder is concerned, any lattice type 
becomes equivalent to a Bethe lattice: In the latter case 
the topology of the lattice allows for retraceable paths 
only, while in the former the statistical properties of the 
J D, s themselves select the same type of paths. Therefore 
it is not surprising that, as in Bethe lattices, in the part 
of the kernel proportional to the disordered interaction 
(S l (r)S. i (rO) ca ^ = (S l (r)S I (rO) for any lattice type. 



3. Spin self-energy and self- consistent procedure 

In the large-z limit and in the paramagnetic (spin- and 
replica-symmetric) phase the local effective action Q is 
cast in the simpler form 



A =-2 



1 

1 dr dr'S(r)S(rO/C(r - t') 

2 Jo Jo 



(14) 



The standard EDMFT procedure is to solve this impu- 
rity problem with a kernel /C self-consistently depend- 
ing on the spin S = Sg via the local susceptibility 
Xioc = Xu an d an (approximate) evaluation of Xij enter- 
ing Eq. (fTB"]) . Here we emphasize again that, while the 
functional dependence of K_A.p(u} n ) from the spin suscep- 
tibilities (cf. Eq. IT3|) ) depends on the dimensionality of 
the AF spin fluctuations, the disorder part always yields 

To solve the impurity problem of Eq. (|14|) is usu- 
ally the difficult step of the EDMFT procedure. In the 
present paper we adopt a Quantum Monte Carlo (QMC) 
technique (see subsection II. B) to obtain the local suscep- 
tibility Xiod^n)- The next step is then to suitably mod- 
ify the kernel of Eq. JTJJ to iterate the self-consistency 
procedure. To this purpose we customarily introduce a 
spin- fluctuation self-energy IIq 1 (u> n ), which, owing to the 
large coordination of the original lattice, is taken to be 
momentum independent. Then in momentum space the 
non-local susceptibility Xij can be written as (see Ap- 
pendix B) 



X(q,w w ) = [n 1 K) + J(q)] 



(15) 
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The calculations reported in Appendix C of Ref. [ll| can 
easily be generalized (see Appendix B) to the present 
disordered case, to obtain the self-energy IIq in terms 
of the local susceptibility and the Fourier-transformed 
kernel 



II 1 (u> n ) = K, AF (uJ n ) + 



1 



Xloc(Un) 

The relation for the averaged local susceptibility 

1 



q 



\ " 



^n - 1 K) + j( q ) 



(16) 



(17) 



provides another equation allowing to relate completely 
the calculated local susceptibility and the bosonic bath, 
allowing to determine from xioc(u) a new IC(u) n ). 
The procedure is then iterated until self-consistency is 
reached. 



B. The quantum Monte Carlo approach 

Within our EDMFT scheme, we are faced with the 
problem of solving the model in Eq. (|14[) with an im- 
purity quantum spin coupled to itself via a retarded in- 
teraction kernel representing the surrounding medium of 
the other spins. As announced above, we here adopt 
a numerical QMC technique, starting from a Hubbard- 
Stratonovich transformation of the local partition func- 
tion 



Woe 



= J DrjTrT exp J drrjij) ■ S(r) 



x exp 



rP 



Jo 



dTdT'K,- X {T,T')r)(T) ■ T](t') 



(18) 



where T is the imaginary time ordering and Tr denotes 
a path integral, taken with respect to the spin degrees of 
freedom S(r). This represents the partition function of 
a spin coupled to an effective imaginary-time-dependent 
random magnetic field with a Gaussian distribution. No- 
tice that within the functional-integral formalism, the 
Hubbard-Stratonovich field r\ is a bosonic field. This 
relates previous analyses of Bose-Kondo model s 20 ! 21 to 
our EDMFT analysis of the magnetic model, where the 
spin degrees of freedom of surrounding medium are rep- 
resented by a simpler bosonic bath. 

Following Ref 1 lOt to implement the QMC algorithm, 
we discretize the imaginary-time axis into L small time 
slices so that the time-ordered exponential in the trace of 
Eq. (|18p can be written as the product of L 2 x 2 matri- 
ces using Trotter's formula. In our calculations we take 
the inverse temperature Jo/3 < 120 and L < 256 while 
keeping Jo At — Jq(3/L < 0.25, where Jo is a natural 
magnetic energy scale: In the 3D case we identify this 
scale with J x (see Sect. IV. A), while in the 2D case we 
use the notation J (see Sect. IV. B). 



III. TRANSITION CRITERIA AND PHASE 
DIAGRAMS 

A. AF instability criterium 

The establishment of an AF long-range order through 
a second-order transition is naturally described within 
the EDMFT formalism via a divergent momentum- 
dependent spin susceptibility 



= [II -V = 0;T N ) + J(Qaf)]" 



30 



(19) 



for some specific AF wavevector and a Neel temperature 
T/ v 11 i 22 . This would occur when the Fourier transformed 
magnetic coupling reaches an extremal value, which, e.g., 
on a hypercubic lattice in D = z/2 dimensions, yields 

J (Qaf) = 1-% E?=i cos(Qaf ) = -Jy/z. 

However, the cumulant expansion and the related 
rescaling of J^ = Jij/^fz, discussed in the previous Sec- 
tion and needed to keep the dynamics of the quantum 
spins, would locate the energy of this ground state very 
far (Eq ~ —yfzJ) from the other states forming the 
"body" of the spectrum. Therefore, having in mind that 
after all real systems occur only up to three dimensions, 
one is naturally led to use a bounded band (with semiel- 
liptic and rectangular densities of states to mimic the 
three-dimensional and the two-dimensional cases respec- 
tively). In this case the divergence of the spin suscepti- 
bility is no longer related to some specific wavevector cor- 
responding to an easily recognizable space arrangement 
of the spins, but is rather reached when the border of 
the band reaches some specific value. Here, after a suit- 
able rescaling of the energy units, we set this condition 
to — J(Qaf) = — J = £min (all the z factors have been 
eliminated via the rescaling). Therefore, the condition 
for AF order reads 



iLo 1 (lj = 0;T n ) = j 



B. Spin-glass instability criterium 



(20) 



On the other hand, one can show that the condition 
for the instability of the paramagnetic phase toward the 
formation of a SG phase is given by 

1 = 4EA"n = Q) = JlY. r l * ; ,2 

V V [no x (o) + J(q)] 



' D 



d e ^ 

-j [n - 1 (o,T SG ) + e 



_ _ T 2 dxiocjO) 

2 ~ D dn^ 



(21) 



where the spin-excitation density of states p(e) = 5(e ~ 
J(q)) has been introduced to calculate the momentum 
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sum. In the particular case when the ordered AF cou- 
pling is absent the action is made local by averaging over 
the disorder irrespective of the EDMFT expansion [cf. 
Section II and Eq. (6) in Ref. H|| and the condition (|2Tj) 
reduces to 1 = Jf)Xf oc = JdXIo<t^ This condition was 
established long ago^ for the determination of the SG 
instability for a quantum Sherrington-Kirkpatrick model. 



C. Dimensional analysis 

We now aim to establish which one of conditions (|20|) 
and (|2ip is realized first depending on the dimensionality 
D. This latter determines the specific form of the bare 
density of states of the spins, which, near the (lower) 
band edge e m in = — J behaves like 



p(e) oc (J + e)~ 



(22) 



To proceed we assume that the condition (|20|) is real- 
ized first and we replace 11^ 1 by J inside Eq. (f2"T]) 



1 = J 2 



j 

de— 



jy 



(23) 



The convergence of the integral is determined by the 
band-edge behavior at (e w —J). For D < 4 the in- 
tegral diverges indicating that, for any finite value of 
Jd, the condition Eq. (f2~lj) is satisfied before T is lowered 
to its Neel value when hLT 1 (Tjv) = J. In this case the 
system always becomes a SG before forming any static 
AF order. We expect that, at least for small Jo's, the 
AF order and a replica-symmetry breaking order coexist. 
It should also be considered that first-order transitions 
might occur driving the system AF or SG ordered with- 
out passing through the instabilities determined by the 
above conditions^. To investigate and settle all these 
issues, however, an EDMFT analysis in the presence of 
finite AF and/or SG order parameters should be car- 
ried out. This is technically challenging and beyond the 
scope of our work. On the other hand, for D > 4, the 
integral in Eq. (|23p converges. In this case a minimum 
value of Jd has to be present, above which a SG order 
is established before an AF order parameter is formed. 
Therefore, while for D < 4 even a small amount of dis- 
order is enough to produce a SG phase, for D > 4 the 
pure AF order is stable and a sufficiently large amount 
of disorder is required to break the replica symmetry. 2 - 

Fig. 1 schematically depicts the phase diagrams for 
our model in two, three, and above four dimensions. 
In two dimensions [Fig. 1(a)], the AF long-range order 
never occurs 2 ^ and AF correlations only manifest them- 
selves with the formation of a spin- liquid phase with a re- 
duced spin susceptibility (with respect to the Curie form 
valid at temperatures above J) and with critical spin 
dynamics. This behavior arises below a crossover tem- 
perature Tsl, which has been numerically determined 
via the QMC analysis reported in Section V. This will 




FIG. 1: Schematic phase diagrams of the disordered AF quan- 
tum Heisenberg model in (a) D = 2, (b) D = 3, and (c) 
D > 4. The dashed line Tsl is a crossover temperature below 
which spin-liquid correlation appear. Tsa and Taf are the 
spin-glass and the AF transition temperatures respectively. 
The dot schematically indicates the position at which the SG 
temperature becomes larger than the AF one. The dotted 
line indicates the "causal instability" line occurring in 2D (see 
text). 



also allow for the quantitative determination of the SG 
transition temperature T$g and of the spin-liquid in- 
stability temperature T*. This latter marks the tem- 
perature below which the spin self-energy lip 1 (a;) as a 
function of real frequencies has a positive imaginary part 
Hq 1 ^)" > thereby violating the causality condition 

wiioV)" <o^£ 

In three dimensions Fig. 1(b) shows that, as soon as 
a disordered three-dimensional magnetic coupling Jd is 
turned on, the system undergoes a second-order transi- 
tion to a SG state. This occurs below a transition tem- 
perature T$g, which, according to our criterion above, 
is always larger than the Neel temperature Tjv (in par- 
ticular it is also larger than the AF temperature in the 
absence of disorder In the opposite limit of Jd ^> J 
one naturally recovers the transition temperature for the 
"pure" quantum Sherrington-Kirkpatrick model— i 23 i 29 
and the Tsg/J curve smoothly approaches the straight 
line T SG = 0.142J D . 

We finally notice that the qualitative picture reported 
in Fig. 1(c) for the case above four dimensions, is in 
qualitative agreement with the mean-field calculations of 
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Refs. 3(J, where the competition between a SG and a 
ferromagnetic phase was studied. 



IV. STATIC PROPERTIES 

A. The three-dimensional case 

When the AF coupling is assumed to act in three di- 
mensions, the spin-wave spectral density p(e) takes a 
semielliptic form 



P3c(e) = — n^~J^ 



(24) 



In this case the integral in Eq. IT7|) can be analytically 
performed to show that this equation is identically satis- 
fied for 



with 



n 1 ( CJ n) = JxXloc{u n ) 



J x = J 2 /4 + Jf 



(25) 



(26) 



Therefore in three dimensions the problem with mixed 
AF and disordered coupling becomes formally equivalent 
to a model with a rescaled coupling J x only 31 . 



1. Phase diagram 

Starting from the free-spin high-temperature regime, 
upon decreasing T, the magnetic correlations start to re- 
duce the magnetic susceptibility. We therefore introduce 
a crossover temperature T x , which is arbitrarily defined 
as the temperature at which the local magnetic suscepti- 
bility is reduced by ten per cent with respect to the free 
moment Curie law 4T X x; oc (T x ) =0.9. From the numer- 
ical calculations, we find T x w 0.75J X . Since we decided 
to use the AF magnetic coupling J as the energy unit, 
we recast this relation as 



(27) 



where the disordered part of the magnetic interaction is 
made explicit. 



2. Spin-glass temperature 

Using the three-dimensional density of states Eq. (|24|) , 
the SG criterion Eq. JH]) is 



2 J 2 



1 = -- 



*J 2 J-j [n - V = 0,T SG ) + e] 



(28) 



Using the relation 



dx = 



the SG criterion can be written as 



v - xn/ 2 



(29) 




-1 = 2 



D 
J 2 



n, 



-i 







J 




Finally, this criterion can be cast into the form 



ru 



J 



1 + 2(J D /J) 2 
y/l + 4(J D /J) 2 



(30) 



(31) 



It is worth emphasizing again that the self-consistent 
value of n^" 1 (Tsg) to be used in the above relations is the 
same as in the non-disordered system. Therefore, to de- 
termine the 3D phase diagram as a function of tempera- 
ture and disordered coupling Jjj , we self-consistently cal- 
culated with QMC the spin self-energy at different tem- 
peratures for a unit value of J x = J = 1. The results 
are reported in the inset of Fig. 2. Then, for each value 
of Jd /J, the SG instability temperature is obtained by 
equating Uq 1 (T = T S g) with the r.h.s. ofEq.©. While 
Fig. 1 is a sketch comparing the phase diagrams of our 
model in various dimensions, Fig. 2 shows the phase dia- 
gram in three dimensions as numerically determined from 
the procedure outlined above. 

It is worth noticing that in the three-dimensional case 
the regime with strong spin correlations reducing the lo- 
cal susceptibility is severely reduced by the occurrence of 
the SG phase. This is to be contrasted with the substan- 
tial spin-liquid regime found in two dimensions (see next 
Section) 



3. Local-moment formation 

The temperature dependence of the static local sus- 
ceptibility is shown in the inset of Fig. 3 for various val- 
ues of the disorder Jd/J at a fixed value of J x . Again, 
by exploiting the fact that in the three-dimensional case 
the disordered coupling can be "hidden" inside the ef- 
fective coupling J x , we first worked in units of J x , and 
calculated once the universal curve JxXioc(T, Jd, J) — 
f(T/J x j). Then we expressed J x in terms of Jd/J 
to find the different curves at various values of disor- 
der in the (meta)stable paramagnetic region, both above 
and below Tsg- It is quite evident that the suscepti- 
bility curves display both a high-temperature and a low- 
temperature Curie- like behavior with a constant (unity in 
log-log plot) universal slope. This is an indication of local 
moment formation, as shown in the main frame of Fig. 3. 
Here the temperature dependence of the square local 
moments H 2 {T) = Txioc(T) is displayed, showing that 
/i(T) 2 w 3/4 at high temperature, while p{T) 2 w 1/4 
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FIG. 2: Phase diagram of the DQAFH model in three dimen- 
sions as obtained from the QMC calculations at fixed (unit) 
value of J x . The solid line represents the SG transition tem- 
perature Tsg as numerically obtained from the criterium of 
Eqs. (f2Tj) and Eq. (f3Tjl . The dashed line is the crossover 
T x (see text). Inset: zero-disorder spin self-energy FI,^ 1 as 
a function of temperature. The antiferromagnetic instability 



occurs at T ~ 0.28J, when II 



o 



J. 



0.75 




0.25 



0.001 



FIG. 3: Temperature dependence of the square effective 
moment /i 2 (T) = Txioc{T) in 3D at fixed (unit) value of 
J x and various values of disorder: Jd/J = 0, 0.5, 1., 5, 10 
from left to right. The high-temperature limiting value is 
S(S + 1) = 3/4, while the low-temperature limit is 1/4. The 
diamonds show the square moment at the corresponding Tsg- 
Inset: Temperature dependence of the local spin susceptibility 
for Jd/J = 0, 0.5, 1, 5, 10 from top to bottom. 



at low temperature. It can also be noticed that the dis- 
ordered coupling strongly favors the formation of static 
low-T local moments. However, as shown by the dia- 
monds in Fig. 3, the SG transition already occurs at 
much higher temperatures when the local moments are 
still close to their high-temperature value 3/4. 



B. The two-dimensional case 

We also consider here the case of two-dimensional 
AF fluctuations having a spectral density of the form 
P2c(e) = jj for — J < e < J- Contrary to the above 3D 
case, the effect of disorder can no longer be included in 
an effective coupling J x together with the AF J. There- 
fore, for any value of Jd one has to recalculate the self- 
consistent quantities xioc, JC- These quantities are re- 
lated by the selfconsistency Eq.(fl"7|). which in two di- 
mensions, reads 



II7V,T) = J coth [Jxioc(u, T)] 



Phase diagram 



(32) 




The results of our EDMFT calculations are summa- 
rized in the numerical phase diagram of Fig. 5. Here 
the presence of various regions is apparent, where the 
system displays different physical properties. First of all 
one should notice a crossover line separating the free-spin 
high temperature region from a region at intermediate 
temperatures and low disorder, where the spin system 
has non-trivial correlations. This crossover temperature 
has been determined from a high-temperature expansion 
of Eq. (j3"2"|) , where the explicit form of n^" 1 in 2D has 
been inserted. Specifically, by assuming that Jxioc is 
small at high temperatures, one obtains that 



(33) 



The prefactor A is then numerically determined, by im- 
posing that the spin susceptibility is reduced by ten 
percent from its high-energy Curie law. This yields 
A w 0.25, for which one obtains the dashed line in Fig. 
4. Below this line and before the SG phase takes place 
at low temperature, the spins form a correlated "liquid" 
with the static and dynamic local susceptibilities show- 
ing anomalous power-law behavior. This behavior, at 
intermediate temperatures where the effects of Jd are 
not significant, is quite similar to the case of the pure 
two-dimensional system described in Ref. [l~3l . Specifi- 
cally, one finds over an order of magnitude in T, where 
the static local spin susceptibility displays the quantum- 
critical power- law behavior Xloc(T) J-^J/T) 2 / 3 . Re- 
lated critical dynamical properties are found in the same 
region (see below). However, upon lowering the temper- 
ature, the paramagnetic phase is characterized by the 
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0.001 0.01 0.1 1 10 

JJJ 

FIG. 4: Phase diagram of the DQAFH model in two dimen- 
sions as obtained from the QMC calculations. The dashed 
line represents the analytic estimation of the crossover tem- 
perature T x the high-temperature free-spin regime to the 
intermediate-temperature spin-liquid regime as obtained from 
the high-temperature expansion (see text). The filled dia- 
monds represent the SG transition temperature Tsg as nu- 
merically obtained from the criterium of Eq. (|2ip . The heavy 
dashed line represents the asymptotic (Jd > J transition line 
for the quantum Sherrington-Kirkpatrick model. The dot- 
dashed line corresponds to the temperature T* of the causal- 
ity violation occurring at zero disorder, while the thin solid 
line is obtained from the low-temperature estimation of Tsg- 

gradual formation of local moments, whose temperature 
dependence for various values of the disorder is reported 
in Fig. 5. 

2. Local-moment formation 




T/J 

FIG. 5: Temperature dependence of the square effective mo- 
ment /x 2 (T) = Txioc(T) in 2D at various values of disorder: 
Jd/J = 0, 0.10, 0.58, 0.98, 2.06 from left to right. The high- 
temperature limiting value is S(S + 1) = 3/4, while the low- 
temperature limit is 1/4. The circle in the Jd/J = curve 
shows the square moment at T = T* « 0.03J. The diamonds 
show the square moment at the corresponding Tsg- Inset: 
Temperature dependence of the local spin susceptibility for 
Jd/J = 0, 0.58, 0.98, 2.06 from top to bottom. The dashed 
line is a guide for the eyes corresponding to a power-law be- 
havior with an exponent 2/3 to illustrate the intermediate 
spin-liquid behavior. 



SG instability as determined by Eq. (|2ip . In the present 
two-dimensional case, the relation Eq. (|2 1 1) characteriz- 
ing the SG temperature can be simplified to 



The fact that disorder favors the formation of local 
moments at progressively larger temperatures, is made 
apparent by the flattening of the curves to the classical 
value 1/4. It is worth noticing that the tendency to form 
a finite local moment at low temperature is also present 
in the absence of disorder (Jd = 0). This tendency was 
also just visible in our previous worki 3 -, in our QMC cal- 
culations at T = J/100. However, it was impossible to 
confirm this numerical hint by well-converged QMC cal- 
culations at lower temperatures 3 - 2 -. On the other hand, 
the static moment formation becomes quite clear in the 
presence of disorder. 



no 1 (w = Q, T SG ) 



J 2 



J 2 = J 



(34) 



with II (ui,T) given by Eq. I|32p In the weak disorder 
limit Jd << J, we approximate xioc by its zero disorder 
expression Jxiociv = 0,T « J) « 0.338( J/T) 2 / 3 (see 
Fig. 1(b) of Ref. 1130 . This gives the asymptotic behavior 
in the limit Jd/J — » 



0.338 

Jd 



3/2 



(35) 



3. Spin-glass temperature 

So far we have discussed the properties of the param- 
agnetic phase disregarding the possible occurrence of a 



showing that a significant Tsg can be obtained even for 
a small disorder. The numerical data seem indeed to be 
smoothly extended to this low-disorder behavior (see Fig. 
4). 
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FIG. 6: Imaginary-time dependence of the local susceptibility 
in three dimensions at different values of the inverse temper- 
ature /3J X = 1, 3, 5, 10, 30, 60, 120. 



V. DYNAMICAL PROPERTIES 
A. Dynamical properties in 3D 

From our QMC analysis we also obtained dynami- 
cal quantities in imaginary time. Fig. [6] displays the 
imaginary-time dependence of the local spin suscepti- 
bility XIoc{t) in the paramagnetic phaseM> Indeed, al- 
though the system rapidly passes from a paramagnetic 
high temperature phase to a low-temperature SG phase 
(see Fig. 2), we follow the paramagnetic phase down to 
low temperatures to allow a comparison with the two- 
dimensional case, where the critical SG temperature is 
much lower. In this way we can investigate the role of 
dimensionality in determining the spin dynamics. 

I n XIoc(t) one can recognize the different temperature 
regimes found in the static quantities and explored in de- 
tail in the previous study of the purely disordered case in 
Ref. This is quite obvious since, as already discussed 
for the static properties, our system in three dimensions 
is equivalent to a purely disordered Heisenberg model 
with an effective magnetic coupling [see Eg. ([25])]. This 
equivalence turns out to be quite important, since the 
dynamical analysis involves a delicate procedure of ex- 
tracting real-frequency information from data obtained in 
imaginary time or frequencies. This is usually a difficult 
task, but it is greatly simplified by the knowledge of well- 
established limits at low and high temperatures for the 
dissipative part of the local response xi oc (lo). This allows 
the construction of an interpolating function that suc- 
cessfully matches the numerical QMC data. Specifically 
one finds a high-temperature regime, where the spin dy- 
namic absorption produces a Gaussian quasi-elastic peak 



xlM _ ps(s + i) 



1-S(S + 1) 



(pj* 



18 



(36) 



with a single characteristic relaxation frequency lot given 
by u\ = 2J X 2 S{S + 1)/2[1 - 5(5 + 1)(/3J X ) 2 /18] In 
imaginary-time this corresponds to 



3 



S(S + 1) 



*l(t) 



(37) 



On the other hand, at low temperatures the curves in 
Fig. 6 clearly display two regimes: At long times the 
curve flattens indicating the formation of nearly static 
moments, while at short times a rapid variation (de- 
crease) of XIoc{t) indicates that some spectral weight is 
shifted at high frequencies. From an analysis similar to 
the one of Ref. QjJ one can see that this leads to dissipa- 
tion from transverse spin excitations 



TTLO 2 



P 



2lot 



cxp 



P(u> — LOtY 



3/2 wtanh(/3w/2) 
V2¥(l + /W2) 
P{lo + lo t ) 2 



Aloi 



exp 



(38) 



at frequencies lot ~ Jy 2 /T. In imaginary time this gives 
T S{S+l) 1 + f3iu T / 2(1 -2t//3) 2 



Xloc( T ) = 



1 + l3u T /2 



x exp 



Plot 



1-1- 



2r 
J 



5(5 + 1) 



$ T (r) 



(39) 



The longitudinal fluctuations instead still produce a 
quasi-elastic absorption peak of the form of Eq. ([3"6"| . but 
with a reduced frequency u>l oc T and a prefactor /35 2 /3 
replacing PS(S + l)/3. This reduction of the effective 
Curie constant is natural since the formation of static 
moments must coexist with the high-energy absorption. 

The known behavior in the high- and low-temperature 
regimes allows to construct an interpolation function 



Xlo 



5(5 + 1) 



[p$i(r) + (1 -p)$ T (r)] 



(40) 



Obviously the temperature variation of the relaxation 
function F(lo) = x' '( w ) / ' {^XiodT)^) is quite similar to 
the one shown in Fig. 2 of Ref. [jjj provided that J is 
replaced by J x • 



B. Dynamical properties in 2D 

Fig. 7 displays the results of our QMC calculations for 
the imaginary-time dependent spin susceptibility. The 
frames from top to bottom are for different values of dis- 
order Jo — 0.1J, Jd — J, and Jo = 10 J. The various 
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curves are calculated at different values of the inverse 
temperature (3J — 1, 3, 5, 10, 30, 60. Frames (b) and 
(c) also contain the curves for (3J — 120, which could 
not be obtained for the zero disorder case due to non- 
convergence of the QMC code. Similarly to the 3D case, 
the low-temperature curves clearly display a flattening at 
large time marking the tendency to form static local mo- 
ments. This formation is clearly favored by disorder. The 
same conclusion could qualitatively be drawn from the lo- 
cal susceptibility as a function of Matsubara imaginary 
frequencies for various values of the disorder: At low tem- 
perature the more disordered systems display a visible 
tendency to form an additional contribution at zero fre- 
quency, which cannot be obtained from a smooth extrap- 
olation from the finite-frequency data. This is because 
the local moments provide a substantial delta-like contri- 
bution to the zero-frequency susceptibility (however, due 
to the difficulty in reliably extracting this additional con- 
tribution, we refrain from showing curves for Xloc(u n ))- 
For the curves of \ioc as a function of imaginary-time, 
we performed an analysis similar to the one carried out 
in Subsection V.A for the three-dimensional case. 

Upon inspecting the two-dimensional phase diagram 
of Fig. [5] one sees that a substantial region is present 
between T x and T* or Tsg, where strong spin corre- 
lations should give rise to a non-trivial spin dynam- 
ics. From our previous work on the non-disordered 
Heisenberg model in two dimensions^, we learned that 
a spin-liquid phase arises below T x . At long imaginary 
times this phase has a critical susceptibility of the form 
XIoc(t) ~ [szn(7TT//3] -1 / 3 , which corresponds to a scaling 
of the imaginary local susceptibility^ 3 . 

XtocM ~ J-Hu/J)- s F s {u/T) , (41) 

with T s {x) = x s \T(^ + i|-)| 2 s inh(f) and 6 = 2/3. 
Eq. (|41[) gives a power-law divergence at low energy, 

X 'U T = M ~ J-*w-§ and xL(?> = 0) ~ 
J-^T-i. 

The question then naturally arises, whether or not this 
critical behavior can be experimentally observed in mag- 
netic systems with strong quantum fluctuations, weak 
disorder, and strong anisotropy in the AF fluctuations. 

In principle one could answer this question by extend- 
ing the analysis of Subsect. V.A and generalizing Eq. 
(|40p to a form including an additional spin-liquid contri- 
bution 

Xioc = — - \p$l(t) + g$st(r)(l -p- «)*r(r)] ■ 

(42) 

with §sl(t) oc [sin(iTT / (3] 1 / 3 . Starting from the 
Fourier transform of Eq. (|4"2")l . we tried to use a standard 
procedure based on Pade approximants to analytically 
continue the Matsubara susceptibility to real frequencies. 
Unfortunately, the analytic continuation procedure is so 
delicate that reliable results can hardly be extracted from 
a form like Eq. (|4"2"|) . In particular we noticed that small 
variations in the relative weights p, q as well as in the 
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FIG. 7: Imaginary-time dependence of the local susceptibil- 
ity in D — 2 for different values of the inverse temperature 
[3J = 1, 3, 5, 10, 30, 60, 120 and three values of disorder: (a) 
Jo = 0.1J, (b) J D = J, and (c) J D = 10 J. For the low- 
disorder case (a) the lowest-temperature data j3J = 120 are 
not available (see text). 



other physical parameters ujl and lot lead to substantial 
variations in the real-frequency susceptibility. Therefore, 
while our analysis is fully compatible with the possibility 
of intermediate frequency windows with anomalous spin- 
liquid exponents, we cannot provide reliable bounds to 
the observation of these regimes. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we systematically investigated the in- 
terplay between disorder and fluctuations in a quantum 
spin-i Heisenberg model. Without entering any broken- 
symmetry phase and neglecting the possible occurrence 
of first-order transitions, we investigated the properties 
and the instabilities of the paramagnetic phase. We 
found that a marked difference exists between the two- 
and the threedimensional cases. In 3D the system is or- 
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dered at a finite Neel temperature Taf for zero disorder 
and becomes a SG at T$g > Taf as soon as a finite dis- 
order is introduced in the magnetic coupling. This result 
is not so surprising because the AF and the SG phase are 
not competing so that antiferromagnetism can well coex- 
ist with glassiness at least for small disorder. According 
to these indications, in real AF systems the presence of 
a small amount of disorder is expected to lead to weak 
glassy properties. These, however, can likely be masked 
by the dominant (at weak disorder) AF properties. 

The two-dimensional case is richer. In the absence of 
disorder, the spin liquid phase analyzed in Ref. [TH was 
shown to be unstable below a temperature T* ~ 0.03 J, 
where the causality condition ujHq (u>)" < on the self- 
energy of the dynamical local spin susceptibility starts 
to be violated. As discussed in detail in Ref. |3J for the 
t — J model, this behavior occurs when long ranged space 
correlations arise at finite temperature being driven by 
a divergent magnetic coherence length £. In this case 
in 2D the real part of the local susceptibility diverges 
because the momentum-independent spin self-energy in- 
troduced in the EDMFT approach is unable to pro- 
duce anomalous dimensions (and the related r\ exponent). 
Then xj oc (0) ~ log£ (unphysically above the bare Curie 
value 35 S(S + 1)/3T). This growth at low (but still finite) 
temperatures entails too large values of the imaginary 
part of x'i'ocj which give an unphysical sign of the spin self- 
energy determined by the self-consistency condition Eq. 
(|3"2")l . This occurs at the "causal instability" temperature 
T* . Nevertheless, the system may display an anoma- 
lous behavior of the spin susceptibility over broad fre- 
quency and temperature ranges before the instability oc- 
curs. Therefore the EDMFT treatment keeps a valuable 
physical interest to detect non-trivial spin correlations in 
low dimensions. Furthermore, we showed here that the 
causal instability may be prevented by a small amount of 
disorder: As seen in Fig. 4, already at Jjj « (0.05 — 0.1) J 
the SG transition temperature overcomes T*. This also 
indicates that EDMFT acquires further interest and va- 
lidity to treat those models where physical mechanisms 
prevent the large growth of the local susceptibility or in- 
duce the transition to a stable phased. 

A noticeable finding of this paper is also related to the 
specific form of the EDMFT kernel obtained in Eq. (fT3)l . 
In particular, as shown in Subsect. II. A and in Appendix 
B, the kernel separates into a AF and a disordered part 
(although disorder is effectively present in both terms via 
the averaged Xij ' s ) an d moreover the disordered kernel is 
build up with the local susceptibility only 

JCd = (Jd/z)J2xu ■ (43) 

i 

The crucial ingredient leading to these two results were 
the local nature of the disorder ((JB Jki)dis = fii,i5j,kJo) 
and the large coordination of the lattice (z — > oo). It 
is interesting to notice that the specific form in Eq. ([43|) 
of the disorder-kernel holds for any lattice type and is 
clearly reminiscent of the form obtained in the case of 



an ordered system on a Bethe lattice^. Formally this 
arises because K,£> is only non- vanishing for (jP jR^dis — 
8i,i5j,kJi), which also holds in the specific case of the 
Bethe lattice where only fully retraceable paths are al- 
lowed. One consequence of the generic validity of Eq. 
(j43| is that the results of Bray and Moore in Ref. [23| 
are naturally recovered here (both in two and three di- 
mensions) in the limit of zero AF coupling. Indeed the 
(cubic) lattice with infinite-range interaction of Ref. 
becomes equivalent to a lattice with infinite coordination 
because the role of the thermodynamic limit (N — > oo, 
with TV number of sites) is replaced here by the z — > oo 
limit. Since the specific form of the lattice becomes im- 
material for Gaussian (local) disorder, our results in the 
Jd/J — > oo limit where K. rs K-d recover both the varia- 
tional and the QMC results of Refs. [I0II23I . In particular, 
as shown in Sect. IV, we find the known linear dependence 
Tsg — 0.142Jd both in 2D and 3D. Our results, however, 
also show that in the presence of a uniform AF coupling 
the SG transition temperature occurs at values larger 
than the above linear behavior. Physically this occurs 
because quantum and thermal fluctuations are naturally 
decreased and severely suppressed by the occurrence of 
the AF phase in 3D. In other words, once AF correlations 
nearly quenches the spins, a (even small) disorder Jo eas- 
ily produces a glassy state at finite temperature Tsg I n 
2D this mechanism is still present but weaker, because 
the AF order does not establish and the spin fluctuations 
are only reduced by the spin-liquid correlations. Then 
Tsg vanishes for Jd tending to zero, but slowly because 
of the logarithmic dependence in Eq. (f35|) . 
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APPENDIX A: DERIVATION OF THE LOCAL 
ACTION 

1. The replica trick 

We consider the DQAFH model defined by the Hamil- 
tonian Eq. {T), and we consider the effect of disorder 
within the replica trick^i. The replica trick relies on the 
assumption that the free energy T of the system is a 
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self-averaging quantity. This means that in the thermo- 
dynamic limit the free energy does not depend on the 
configuration chosen for the random couplings jP. As a 
consequence, T remains unchanged when it is averaged 
with respect to all the possible configurations of J 9 , and 
we find 



(3T 



(logZ(jP) di 



(Z{JE))di. 



- lim - — 

= lim hog(Z n (jP)) 
n^o n 



dis 



(Al) 



where Z(J®) is the partition function of the system for a 
given configuration of J® , and (• • • )dis denotes an aver- 
age with respect to the probability distribution P(J®). 

One of the replica trick hypotheses consists of assuming 
that the average (Z n (J®))di S is an analytical function of 
n. We can thus express this quantity for integer values 
of n, interpreting Z n (J®) as being the partition function 
of n replicas of the system which are independent of each 
other but characterized by the same distribution of J®. 
This is formally expressed by the relation 



{z n {JS))* 



n 



+oc 



dJ»P{J»)Tr 



exp - 



(A2) 

where H a is the disordered Heisenberg Hamiltonian 
Eq. (JXJ) acting on the replica a. In Eq.([3]), the trace Tr is 
applied to all the operators Sf of the replicated Hilbcrt 
space. Here, the integration with respect to J? gener- 
ates an effective action with the following properties: 

(i) The corresponding effective model now recovers the 
translational symmetry of the underlying lattice. This 
allows a dynamical mean field (DMFT) approach. 

(ii) The spins of the same replica are coupled through 
the antiferromagnetic exchange J^ F ', and an antiferro- 
magnetic (AF) instability can occur. 

(Hi) Some correlations occur between different replicas 
of the same spin, which could give rise to a ground state 
characterized by a replica symmetry breaking. This state 
would describe a SG phase. 



with 



V< 



jAF 
J 0j 



J 0j 



(A4) 



describing the magnetic interaction between the site i = 
and its neighbors. The usual static mean field approxi- 
mation consists of replacing this term by its average with 
respect to the thermal and disorder fluctuations. In a 
paramagnetic phase, this contribution would vanish and 
we would have to go beyond the static mean field, within 
the DMFT approach. Using a path integral formalism, 
the partition function of the replicated system is given 
by 



Z n (J°)=TrTexp 



cItH(t) 



(A5) 



where T is the imaginary time ordering and the trace 
Tr denotes a path integral, taken with respect to all the 
spin degrees of freedom Sf (r) . This is now formally per- 
formed in two steps Tr — TroTA°\ by considering first 
all the degrees of freedom concerning the sites j ^ 
(Tr^), and then computing the partial trace for the lo- 
cal spin operators on site i = (2>o). We find 



Tr Tr^ T exp 



(A6) 



We average this relation with respect to the random dis- 
tribution of the couplings JP , assuming that the proba- 
bility distribution P(J®) has no correlation between two 
different pairs of sites \i, j). The averages for the cavity 
(with £f(°)) and the local (with Voj) parts can thus be 
factorized 



(Z n (jP)) dls =Z<- 



x (Tr (T exp 



f"drV 0j (r) 



\(0)' 



dis 



(A7) 



where 



2. DMFT and the local partition function 

Following the usual DMFT approach^, we consider 
the local effective partition function seen from a given 
site of the lattice, that we arbitrarily denote i = 0. As a 
preliminary, we introduce the 'cavity' Hamiltonian , 
corresponding to the Hamiltonian H without the site i = 
0. The complete Hamiltonian is related to the cavity one 
with 



H = 



z (o) _ (7Y(o) Texp 



drH W (r) 



) dis , (A8) 



and (■■■)^ denotes the thermal and disorder average 
with respect to the cavity Hamiltonian This is 

performed by using a cumulant expansion, using the gen- 
eral formula (Texp[A])<® - Texp\Z£S C P /p% where 

A = — J>2j*to Jo drVoj (t) and C p is the p— order cumu- 
lant, which invokes correlation functions between p op- 
(A3) erators S". The usual static mean field approximation 
would consist of neglecting all the cumulants p > 2. But, 
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since here we consider the paramagnetic phase, all the 
cumulants C p with odd p vanish. Furthermore, in the 
limit z^oowe find C p m z 1 ^/ 2 , and (Texp[A])W = 
Texp[C2/2]. The DMFT approach, which is exact in the 
limit of a large coordination number z, thus consists of 
approximating the relation Eq. (|A7|) with 

(Z n (J§)) dis = Z< ><Tr <Texp^>< >U . (A9) 



2. Disorder averaging of the kernel 



where 



Ci=\ dr f dr' E w -EwS^(r)^(r') 
J a Jo 

x £«^<ts£(t)s£,(t')>< > 



(A10) 

From this expression Eq. (fj| is then obtained, which 
reduces to Eq. ([13)1 in the paramagnetic case. 



APPENDIX B: DISORDER AVERAGE OF THE 
ACTION KERNEL 

1. Preliminaries 



Here we will average the large— z expansion (|B2p with 
respect to all the possible distributions of J® . First, we 
remark that a given bond I — m (and the related Ji m ) 
either appears in the product of J's, forming the "skele- 
ton" of the path, or is inside the loops forming one (and 
just one) of the IPs. As a consequence, when taking the 
average over the disorder, each object, a J or a II is sep- 
arately averaged with respect to the others; therefore the 
average over the disorder distribution P(j(j) factorizes in 
the average of each separate II and of the product of J's. 
Furthermore one should notice that, since the paths are 
non-self-intersecting, the product of J's can only contain 
one J per I — m bond. Therefore, owing to the zero av- 
erage of the disordered J D , this product vanishes unless 
it only contains purely AF ordered couplings. Then 

Xij = SijILii ~ -^IluTljj 



E 

paths 



(_-\\p jAF jAF jAF 
\ L ) J ii 1 J i l i 2 ' ' ' J i p j 



(B3) 



Our treatment of disorder proceeds along the following 
steps: First we fix the disorder configuration assuming a 
specific realization of the Jy distribution. Then, for this 
configuration one can closely follow the same steps as in 
Ref. [ill to derive the relation between the spin correlation 
functions with and without cavity 



X 



cavity 

ij — 



Xij — XiOXOj/XOO ■ 



(Bl) 



Here, we want to establish a similar relation for the 
disorder-averaged susceptibilities Xij defined by Eq. (J9j) . 

In general one can see that the various correlation func- 
tions are obtained by sums of the type 



Xij 



^r/Hu . — n^Hy 

V Z 



E 

paths 



( 1) P Jtii Ji±i 2 ' ' ' Jipj 
z (P+l)/2 



(B2) 



where, to obtain the leading contribution in the large- 
z expansion, the sum is taken over non-self-intersecting 
paths and the local correlations n take into account loop 
decorations of these paths. We will see in the following 
how these local propagators can be related to a local self- 
energy, which will be determined in a self-consistent way. 

The relation (|Bip for the cavity susceptibility is ob- 
tained from Eq. (|B2|) by using the general algebraic 

identity Ei^o Eo^j = Ei-j ~ YTHf V ] Eo^o> where 
Ej— >j represents a summation over all the direct paths 

between two sites i and j, and Y^i°^j V excludes all the 
direct paths through the site 0. 



where the bar indicates disorder-averaged quantities. 
Then the same arguments used in Ref. [ll| for obtaining 
the relation (|B1[) from the expansion (|B2[) . can now be 
applied to the averaged quantities of Eq. (|B3j) giving 



— cavity — — — / — 

Xij — Xij ~ XioXoj I Aoo 



(B4) 



This completes the derivation of Eq. ([IT]) leading to the 
final expression of the kernel Eq. (fT3|) 



3. Local self-energy 



The expansion (|F33|) is equivalent to a Bethe-Salpeter 
equation for the average susceptibilities 



AF 



,(w„) = <%n (w„) - n (w„) ^2 -^xiji^n) , (B5) 



where we restored the explicit Matsubara frequency de- 
pendency, and the averaged n = Ha is site independent. 

Averaging the distribution of J® restores the trans- 
lational symmetry of the underlying lattice. We can 
thus define the q— dependent susceptibilities x(Qj w n) as 
the spatial Fourier transform of ^-(w,,). Introducing 
J(q), the Fourier transform of the magnetic exchanges 



z, the Bethe-Salpeter Eq. (|B5j) can be written as 



X(q,Wn) = [n + J(q)]' 



(B6) 



The AF part Kaf of the Kernel is expressed in terms 
of the averaged susceptibilities using the definition [cf. 
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Eq.®] 



Using this relation together with Eq. (|B6[) . and the iden- 
tity Xioc{u n ) = J2qX(l^n), we find 



-Xio( W n)Xoj(^")/XooK0] • (B7) 

We rewrite this expression in momentum space: 

K,AF{u n ) = ^ J 2 (q)x(q,o; n ) 
q 

-i 2 



II (u>„) = K.AF(^n) + 



(B9) 



J (q)x(q, 



/XiocK) . (B8) 



An explicit similar calculation can be found in Appendix 
C of Ref. Gil 
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